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Abstract. The excitation rate P of solar p-modes is computed with a model of 
stochatic excitation which involves constraints on the averaged properties of the solar 
turbulence. These constraints are obtained from a 3D simulation. Resulting values 
for P are found ~ 9 times larger than when the calculation assumes properties of 
turbulent convection which are derived from an ID solar model based on Gough 
(1977) 's formulation of the mixing- length theory (GMLT). This difference is mainly 
due to the assumed values for the mean anisotropy of the velocity field in each 
case. Calculations based on 3D constraints bring the P maximum closer to the 
observational one. 

We also compute P for several models of intermediate mass stars (1 Af < 2 Mq). 
Differences in the values of Pmax between models computed with the classical mixing- 
length theory and GMLT models are found large enough for main sequence stars 
to suggest that measurements of P in this mass range will be able to discriminate 
between different models of turbulent convection. 
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1. Introduction 

Excitation of solar-type oscillations is attributed to turbulent move- 
ments in the outer convective zone of intermediate mass stars. 

Accurate measurement of the rate P at which acoustic energy is 
injected into such oscillations will be possible with the future seismic 
missions (e.g. COROT and EDDINGTON). Comparison between mea- 
sured and theoretical values of P obtained with different models of 
turbulent convection will then provide valuable information about the 
properties of stellar convection zones. 

Models for stochastic excitation have been proposed by several au- 
thors (e.g. Goldreich and Keeley, 1977; Balmforth, 1992; Samadi and 
Goupil, 2001). In the present work, we consider the formulation of 
Samadi and Goupil (2001, Paper I hereafter). Constraints on the time 
averaged properties of the solar turbulent medium are obtained from 
a 3D simulation. They allow us to compute P and compare it with 
solar seismic observations and results obtained with Gough (1977)'s 
formulation of the mixing-length theory (GMLT hereafter). 




© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



samadi.tex; 2/02/2008; 0:51; p.l 



2 R. Samadi 

2. The model of stochastic excitation 



According to Paper I, the rate at which a given mode with frequency 
cuq is excited can be written as: 

, f^. ^ 4 /16 ^4 /a,s\2 5 1 

P(a.o)ocy^ dm-pow [j-^^ [^) ^« + 3 j 

(1) 

In Eq. (1), pq is the mean density, is the radial component of the 
fluid displacement adiabatic eigenfunction ^, as = (dp/ds)^ where p 
denotes the pressure and s the entropy, is the rms value of the 
entropy fluctuations, gr{S,r, fn) involves the first and the second deriva- 
tives of ^ with respect to r, Sji{iOo,m) and Ssic^o^i^) driving 
sources inferred from the Reynolds and the entropy fluctuations respec- 
tively, $ is Gough (1977)'s mean anisotropy factor defined as $(m) = 
<u^ > — < u>'^/w'^ where u is the velocity field, < . > denotes time 
and horizontal average and w is the mean vertical velocity (w"^ =< 
u1 > — < Uz >^). Expressions for Sii{uo,m) , SsicoojiTi) and gr{^r:,fn) 
are given in Paper I. 

The driving sources Sr and Sg include the turbulent kinetic energy 
spectrum E{k,m), the turbulent spectrum of the entropy fluctuations 
Es{k, m) and Xk{^) the frequency-dependent part of the turbulent spec- 
tra which model the correlation time-scale of an eddy with wavenumber 

Xk{'^) is modeled here with a non-gaussian function constraint from 
the 3D simulation. 



3. Numerical constraints and computation in the solar case 

We consider a 3D simulation of the upper part of the solar convective 
zone as obtained by Stein and Nordlund (1998). The simulated domain 
is 3.2 Mm deep and its surface is 6 x 6 Mm^. The grid of mesh points is 
256 X 256 X 163, the total duration 27 mn and the sampling time 30s. 

The simulation data are used to determine the quantities E{k,m), 
Es{k, m), w, s and $ involved in the theoretical expressions for Sr, Ss 
and P. More details will be given in forthcoming papers. 

We compute P according to Eq. (1): the cigcnfunctions (^) and their 
frequencies {u = ojq/2ti) are calculated with Balmforth (1992)'s non- 
adiabatic code for a solar ID model built with the GMLT approach. 
The /c-dependency of E{k,z), is modeled as (k/ko)^^ for k < kg and 
as {k/ko)~^^^ for k > ko where ko = 27r//3A, A = aHp is the mixing- 
length. Hp the pressure scale height. The value of the mixing- length 
parameter, a, is imposed by a solar calibration of the ID GMLT model. 
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The value of ko - hence of /? - is obtained from the 3D simulation. The 
above analytical /c-dependency of E and Es{k, z) reproduce the global 
features of E and Eg derived from the 3D simulation. 

Results of P computations are presented in Fig. 1 (solid curve) 
and compared with values of P (filled dots) derived from solar seismic 
measurement by Chaplin et al. (1998). 

We also compute P in the case when the quantities involved in 
Eq. (1) are all obtained from the GMLT model. In the mixing- length 
approach, the anisotropy factor, is a free parameter. We then assume 
two different values for $ = 1.3745 (dot dashed curves) and $ = 2 
(dashed curves). The value = 1.3745 used in the GMLT model 
provides the best fit between computed solar damping rates and the 
measured ones by Chaplin et al. (1998). 
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Figure 1. Solar p-modes excitation rate Piv) (See text for details). 

Without any adjustement of scaling parameters and instead using 
all the constraints inferred from the 3D simulation, we find a maximum 
of P much larger (~ 9 times larger) than is obtained using a ID 
GMLT solar model with <5 = 1.37. This difference is mainly due to 
the assumed value for ^ in the excitation region. Indeed, analysis of 
the 3D simulation suggests that $ ~ 2 in the excitation region. When 
$ = 2 is used to compute the excitation rate with the GMLT model, 
P comes close to P calculated with the 3D simulation constraints. 

Our result then shows that the values of $ found for the solar GMLT 
model when adjusted on the damping rates is not compatible with the 
actual properties of turbulent medium in the excitation region. An 
improvement could come from a consistent calculation which would 
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assume a depth dependent $(m), as suggested by the simulation, in 
both damping and excitation rates computation. 

Our calculations using the 3D constraints bring Pmax, the P maxi- 
mum, closer to the solar seismic measurements but still under-estimates 
them by a factor ~ 2. More sophisticated assumptions for ko{z) and A 
will likely lead to a better agreement with the observations. 

4. Scanning the HR diagram 

We consider two sets of stellar models previously investigated by Samadi 
et al. (2001): Models in the first set are computed with the classical 
MLT whereas those belonging to the the second set are computed with 
GMLT. Wc find that the maximum of P, Pmaxi scales as 3.2 log(-L/L0 x 
Mq/M) for the first set and scales as 3.51og(L/L0 x Mq/M) for the 
second one (L, Lq, M and Mq have their usual meaning). This result 
suggests that measurements of P in several differents intermediate mass 
stars (1 < M < 2 Mq) will enable one to discriminate between different 
models of turbulent convection. 
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